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Abstract 

(N 

> 

\Q ' To know the space time evolution of a heavy ion reaction 

o ■ 

' is of great interest, especiahy in cases where the measured 

I spectra do not ahow to ascertain the underlying reaction 

mechanism. In recent times it became popular to believe 
that the comparison of Hanbury-Brown Twiss correlation 
functions obtained from classical or semiclassical trans- 
port theories, like Boltzmann Uehling Uhlenbeck (BUU), 
U ■ Quantum Molecular Dynamics (QMD), VENUS, RQMD or 



ARC, with experiments may provide this insight. It is the 
purpose of this article to show that this conjecture encoun- 
ters serious problems. The models which are suited to be 
compared with the experiments at CERN and Brookhaven 
are not able to predict a correlation function. Any agree- 
ment with existing data has to be considered as accidental. 
The models suited for lower energies can in principle pre- 
dict correlation functions. The systematic error may be 
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too large to be of use as far as quantitative conclusions are 
concerned. 



lypeset using REVI^ 



I. INTRODUCTION 



It is a common problem in heavy ion reactions between 25 MeV/N and 200 GeV/N 
that the single particle spectra do not allow to ascertain the underlying reaction 
mechanism. To mention only two examples: At low energies despite of many years of 
efforts the fragmentation of nuclei into many intermediate mass fragments remains 
still a process whose origin is heavily debated. At high energies it turned out to be 
very difficult to rule out a hadronic scenario which may produce the same spectra as 
those proposed ClS db SI gnal for the creation of a quark gluon plasma. 

In such a situation a search for experimental information beyond the single parti- 
cle spectra is obvious. Most valuable would be an information on the spatial structure 
of the reaction. It would allow to calculate key quantities like densities or energy 
densities. This information is, however, hard to obtain. 

The only promising method proposed up to now is based on the interferometry 
of identical particles. The interference of the amplitudes of two indistinguishable 
processes gives rise to a correlation function which in principle allows to extract the 
radius of the emitting source. This approach has been very successfully applied by 
Hanbury-Brown and Twiss [jT| in astronomy to determine the angular radius of stars 
by measuring the spatial correlations between two photons. Later Goldhaber and 
Kopylov and Podgoretsky advanced its application in particle physics by showing 
that measurable momentum space correlations may contain information on the size 
of the emitting source. 

In the ideal case of a large, randomly emitting source of known shape this method 
is indeed very powerful and the experimental results can be directly related to the 
source radius of the emitting object. In particle and heavy ion physics the situation 
is, however, much more difficult. There we encounter quite a number of problems. 
The size of the emitting sources is of the order of the radius of a nucleus and therefore 
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not small as compared to the size of the wave function of the emitted particles. This 
renders some approximations impossible. The signal may be distorted by final state 
interactions between the emitted particles or by to the long range Coulomb force of 
the source. The emission time point cannot be defined unambiguously. The decay 
of resonances into identical particles or correlations between the momenta of the 
emitted particles and the coordinates of the emission point may pretend a wrong size 
of the source. For a discussion of these problems we refer to ref. [H , @. Recently 
is has been discussed that the HBT correlation function for an expanding source, as 
encountered frequently in heavy ion reactions, yields a much more difficult relation 
between the space time structure of the emitting source and the correlation function 
as that for a static source |^. 

Due to these problems the measured correlation function in heavy ion collisions 
cannot be directly related to the parameters of the emitting source even if its form 
were known. In this situation there are two possibilities. Either one assumes the form 
of the source and uses the measured correlation function to fix the source parameters. 
Unfortunately this procedure makes these parameters model dependent. Hence they 
cannot be used for more than a comparison between different experiments and yield 
little information on the actual source properties. Or one tries to describe the reaction 
in its entity. This turns out to be a quite complicated procedure but became very 
popular recently. In this approach one follows the time evolution of the system with 
help of one of the standard transport models. Unfortunately none of them propagates 
(anti) symmetrized wave functions but at most a direct product wave function. Since 
the HBT effect is based on the (anti)symmetrization of the wave function of identical 
particles the transport model themselves cannot predict the correlation function. 
Rather one assumes that each particle "freezes out" at some time point. The freeze 
out time is different for each particle. At high energies it is assumed that the freeze 
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out time is that time at which the particle encounters its last collision with another 
particle of the system. At low energies, where potential interactions are important 
as well, the freeze out time cannot be unambiguously defined. The freeze out times 
as well as the particle momenta and positions form then the input for the subsequent 
calculation of the Hanbury-Brown and Twiss (HBT) correlation function |p which 
is then compared with experiment. 

Agreement is usually interpreted as a sign that the underlying transport model 
gives a realistic space time evolution of the different particles. Since these transport 
codes provide not only the momentum space coordinates of the particles but also 
that of the coordinate space they can then be used to calculate the time evolution of 
key quantities like the energy density or the density. 

The weak point in this procedure is the transition between the transport model 
and the subsequent program which calculates the correlation function. Does the 
transport model provide the correct time evolution of those quantities which are 
essential for the calculation of the correlation function? 

It is the purpose of this article to show that this is hardly the case. In order to 
understand the reason one has to understand in detail the derivation of the different 
transport models from the fundamental quantal equations as well as the derivation 
of the equation which is employed to determine the correlation function. We will 
perform this investigation for the three types of present day simulation programs: 
The Quantum Molecular Dynamics approach (QMD) 0, BUU type models like 
Boltzmann Uhling Uhlenbeck (BUU) llO|- [|ll|], Vlasov Uhling Uhlenbeck (VUU) 
or Landau Vlasov (LV) ||13|| and cascade models. The later class includes also the 
high energy simulation programs like VENUS |T^, RQMD and ARC [|16|. 



This problem is independent of the relativistic or nonrelativistic nature of the 
approach. It also does neither depend on the time between the emissions of the 
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two identical particles nor on the presence of resonances. It is also independent of a 
possible final state interaction between the particles which is therefore omitted. The 
common demand on all transport programs is that an emission time point can be 
defined. Essential is the information the programs provide at that time point. This 
information is quite different for the three types of programs mentioned above and 
hence the systematic errors are specific to each of the different transport models. In 
two cases (QMD and BUU) this procedure implies systematic errors which question 
the usefulness of the approach for its original purpose: The discrimination between 
different reaction mechanisms which yield the same single particle spectra. For the 
high energy simulation programs the correlation function is completely artificial. 

For clarity we limit our formalism to the simplest form possible by assuming that 
we are dealing with two bosons which are simultaneously emitted and can be treated 
nonrelativistically. For this simple case the formalism is very transparent. More 
realistic but also more complicated scenarios may add additional problems but do 
not overcome the problems discussed here. 



We start with the derivation of the correlation function which relates the freeze 
out points with the measurable two body correlation function. We assume that a 
source, which is considered as classical, emits simultaneously two identical bosons 
with momenta pi and p2- The differential two body probability W reads then as 
follows [0 (;i,c= 1): 



a and (3 characterize the emitting source. The (anti) symmetrized production ampli- 
tude Ts is given by 



II. THE CORRELATION FUNCTION 



dpidp2 
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Ts{pi,p2,a,(3) = -^{T{pi,p2,a,P) ±T{p2,pi,a,pj). (2) 
where T{pi,p2, a, /3) is the Fourier transform of the wave function 

Tip„p2,a,P) = J < > . (3) 

Introducing the Wigner density of the two body density matrix p2 = \ip2 >< "^21 

D{xi,Pi,X2,P2) = J n rfW^^'^^' < ^1 -yi/2,X2 -y2/2|^2 >< ^2|^1 +yi/2,f2 + y2/2 > 



j=l,2 



(4) 



we can express the probabihty as a function of the two particle Wigner density 

-j^—r = d Xid X2[D{Xi,Pi,X2,P2) ± D{xi, — - — ,X2, — - — )cos{pi - P2){xi -X2) . 
apidp2 J II 

(5) 

Hence for the calculation of this probability the transport theories have to provide the 
two body Wigner density. Unfortunately most of them do not permit to calculate 
this quantity. Therefore one has introduced an approximation, called smoothness 
assumption (SA) : 

^J\x\-, — ^ — ,a;2, — ^ — ) ^ Z;(xi,pi,X2,P2)- (6) 

We will discuss the limits of its validity which turns out to be crucial in the course of 
the article. Employing the smoothness assumption the two body probability reads 
as 

d^W^^ 



)\ap2 



d^Xxd^X2D{xx,P\, X2,P2){1 ± cos{pi - P2){Xi - X2)). (7) 



This is the standard expression for the two particle probability employed in numerous 
publications to relate the measured cross section with the radius of the emitting 



source. 



Depending on the available information, for actual calculations one may have to 
employ further approximations. For BUU type equations, which propagate the one 
particle Wigner density only, one assumes that the correlations between particles are 
neghgible 

D{xi,Pi,X2,P2) ~ D{xi,pi) ■ D{X2,P2) (8) 

whereas for classical cascade calculations one assumes that the quantal two body 
Wigner density can be replaced by the classical 2 body phase space density Fd 

D{xi,pi,X2,P2) ^ Fa{xi,pi,X2,P2)- (9) 

For a static source without any correlation between the emission point and the mo- 
mentum of the emitted particle the correlation function, the quantity one compares 
with experiment, is independent of the center of center of mass motion of the emitted 
pair and is given by 

•> dpi dp2 

where P — ^'^+p^) jg the center of mass momentum, p = pi — p2 is the relative 
momentum and ^ is the one particle momentum distribution 

Jd'x,D{x„p^). (11) 

In the general case, where correlations are present, C depends on the center of mass 
motion as well. As we will see the correlation function C{p) contains the desired 
information about the spatial properties of the emitting source. 

III. CONSEQUENCES OF THE SMOOTHNESS ASSUMPTION FOR 

CASCADE CALCULATIONS 



Cip) - rd^Z^^r. (10) 



One of the classes of models employed to extract source radii by comparing ex- 
perimental results with model predictions are the so called cascade models. These 
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are classical n-body models which solve the Hamilton equations of a n-body system 
and are presently the only models available to simulate heavy ion reactions at CERN 
and Brookhaven energies. They include VENUS |T^, RQMD [|T^ (in its usually 



employed cascade version) and ARC |T6[ as well as now less frequently employed 



programs for heavy ion reactions at an energy of around 1 GeV/N [|T9|. In these 
models the particles do not interact via potentials but suffer two body collisions if 
they come sufficiently close in coordinate space. In between the collisions the parti- 
cles move on straight lines. In the computer programs they are treated as classical 
particles with a sharp momentum and a sharp position. 

One may ask how classical models can be employed to calculate a correlation 
function which is solely based on the interference of amplitudes and hence a genuine 
quantal effect. For an understanding we have to make a detour. In order to employ 
eq. 5 we have to construct the Wigner density out of the classical two body phase 
space density. This is of course not unique but the approach 

Fcl = 5{Xi- Xa)5{pi - PI3)S{X2 - X^)5{p2 " P^) 

= lim £!:^g-{P-i?(i))2C/4-(X-i?(t))24/Dg-(p-fc(t))2C-(x-r(t))Vl5 
C^oo,D->0 TT^ 

= Dcl{Xi,Pi,X2,P2) (12) 

serves our purpose. The expression in the last line will be considered as Wigner 
density.Here we have used the definitions 

P = pi + P2; K = Pa + Pp 

_ Xl + X2 _ ^ _ Xg+Xp 

X = Xl — X2] r = Xa — Xp. (13) 
The Wigner density Dci{xi, p^+p^ , X2, ^^^^^ ) is obtained by replacing pi and p2 by 
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Please note that this Wigner density does not respect the uncertainty relation. 
Inserting this expression in eqs. 5 and 7 and performing the limit procedure we 
obtain for the two particle correlator without smoothness assumption 



= 5{p- k{t)) ± 5{k{t))cos{2pr). (14) 



r dW dW 
J dpi dp2 



d^P = Sip- k{t)) (15) 



This expression differs from 

dW dW 

only for the case that the relative momentum of the emitted classical particles is zero 
what in practical terms never happens. For all other cases we find 

Cip) = 1. (16) 
Applying the smoothness assumption we find, however, 

/ " ,^ d'^P = 5(p- k(t))(l + cos2pr) (17) 
J dpidp2 

and hence 

C^^ip) ^l + cos2pr. (18) 

Thus we observe that here the smoothness assumption creates correlations out of 
nothing. One faces the somewhat surprising result that the correlation function and 
hence the extracted radii are artificial and are only due to the differences between 
the approximate and the exact formula for the correlation function. Applying the 
correct formula the cascade calculations do not yield any correlation function, as the 
exact result shows. 

The truth of this observation can even easily be verified without any calculation. 
If two particles with a sharp momentum are emitted from two localized sources one 
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can measure the momentum sufficiently precise in order to identify the source from 
which each particle has been emitted provided the both momenta are not identical. 
Thus there are no alternative processes, hence no interference of their amplitudes and 
there is, as a consequence, no HBT correlation function. This has unfortunately the 
consequence that there is presently no microscopic model which may be used for the 
interpretation of the correlation data measured with ultrarelativistic heavy ion beams 
at CERN and AGS. 



The discussion of models which provide sufficient information to construct a cor- 
relation function we begin with the QMD approach because it is the only one which 
allows to calculate the 2 body Wigner density. Hence one can calculate (cq.5) 
without any approximation. One can furthermore introduce the smoothness assump- 



validity of this approximation in the situation of a heavy ion reaction and hence for 
the judgement of the predictive power of the correlation function calculated in the 
framework of the other models. 

The QMD model is a n body theory which simulates heavy ion reactions between 
30 MeV/N and 2 GeV/N on an event by event basis. Each nucleon is represented by 
a coherent state of the form 



Thus the wave function has two time dependent parameters Xa,Pa, L is fixed. As we 
will see this wave function serves as a test wave function for a variational principle. 
Hence it is an input of the calculation and not the result of the solution of the 
Schrodinger equation. It relies heavily on intuition; other test wave functions may 



IV. QUANTUM MOLECULAR DYNAMICS (QMD) 



tion and can calculate then applying eq. 7. This may serve as a test for the 




(19) 



11 



yield a different time evolution of the system. The total n body wave function is 
assumed to be the direct product of n coherent states 

= (f)a{xi,Xa,Pa,t)(j)l3{x2,Xp,p^,t) (20) 

thus antisymmetrization is neglected. The initial values of the parameters are chosen 
in a way that the ensemble of At + Ap nucleons gives a proper density distribution as 
well as a proper momentum distribution of the projectile and target nuclei. The time 
evolution of the system is calculated by means of a generalized variational principle: 
We assume that Pa and contain the essential time dependence of the n-body wave 
function. The Lagrange function C can then be written as a functional of these 
parameters where H is the n - body Hamiltonian . 



. d ^ dpa d ^ dxa d . ^ 
dt dt dpa dt dxa 



(21) 



The time evolution of the parameters is obtained by the requirement that the action 

t2 



Jc[<P,cP*]dt (22) 



h 

is stationary under the allowed variation of the wave function. For the wave function 
of eq. 20 the Lagrange function is given up to a constant by 

C = UpJa - ^ - ^ - IEV{x^,x,)). (23) 

V{xa,Xi3) is the expectation value of the (density dependent) 2 body potential. The 
variation of the Lagrange function gives Euler Lagrange equations for each of the 6 
parameters 

i^a = -V^.J2Vi^a^^p) (24) 

and 

= Pa/m. (25) 
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With these two equations one has reduced the problem of solving an- body 
Schrodinger equation to that of solving 6 n ordinary differential equations. In reality 
V is a parametrization of the real part of the Briickner G- matrix. The imaginary 
part is approximated by the measured cross section. For details we refer to ref. 0. 
Hence in QMD the centroids of the Gaussians in momentum and coordinate space 
are the only quantities which change in time. The form of the wave function around 
the centroids is fixed. This is a consequence of the ansatz ( eq. 19). 

From the test wave function eq.(19) we calculate the Wigner density of a pair of 
particles 



Inserting this Wigner density in eq. 5 one obtains after integration over the pair's 
center of mass momentum 



where t is the (assumed common) freeze out time. This is the probability to find two 
particles with a relative momentum p, which have been emitted from two classical 
sources at a relative distance of f and a relative momentum of k. 



In order to derive the equation for the time evolution of the one-body Wigner 
density of a particle moving in a selfconsistent potential V{x) we start from that for 
the one body density matrix pi = ><ipi\ 



<'^L~{x-f{t)f/L 



(26) 




(27) 



V. ONE BODY TRANSPORT THEORIES 



Pi 




(28) 



Applying to this equation the Wigner transformation for an operator O 
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1 — # — * 

Ow{x,p) = ^/ d'ye'^Kx- ||0|f + | > (29) 
one obtains the differential equation 

(^ + ^V3,)D{x,,p,,t) = I Sp\K,{p,-p\,x,)D{x,,p\,t) (30) 
D being tlie Wigner density of tlie one body density operator and Ki is defined as 

K,{p, - p[, x) = ^j -0^,e-^^^'-''^^''\V{x + y/2) ' V{x - y/2))- (31) 

We liave restored % liere for reasons wliicli will soon become obvious. One can expand 
the integration kernel around Xi 

KM -P'i,xi) = ^M^^^)V{x,)S{p,-p[). (32) 

We see that Ki can be viewed as a series with the expansion coefficient /?.V^^Vp'i. 
Hence in the limit that the expansion can be terminated after the first term the 
Schrodinger equation in its Wigner representation is equivalent to the classical Vlasov 
equation : 

d vt -> -> -> 

+ -V^jL'(fi,Pi,t) = (V^,V^(fi))Vp.,D(fi,pi,t) (33) 
ot m 

The Vlasov equation describes the time evolution of the phase space density of 
particles which move on classical orbits specified by the Hamilton equations ^ = ^ 
and ^ — —'Vg^V. As in QMD V presents the real part of the Briickner G- matrix 
and the imaginary part is added as a cross section. 

There are two approaches to solve the above equation. Either one solves the 
differential equation directly or one creates a swarm of test particles which are subject 
to the Hamilton equations and fulfil the initial condition D{xi, pi, to). One propagates 
this swarm with help of the Hamilton equations until a time t and then constructs 
the Wigner density D{xi,pi,t) by coarse graining. The latter solution method is 
called test particle method and is employed in the BUU, VUU and LV approaches. 
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When calculating the observables, i.e. the expectation values of operators, the 
transition from the first to the second method corresponds to the replacement of an 
analytical integration by a Monte Carlo procedure. Using the swarm of test particles 
the analytical solution 

{0{t)) = I D{x,,p,,t)0{x,,p,)d^x,d% (34) 

is replaced by the corresponding Monte Carlo type integral 

1 ^ 

1=1 

where the fi{t) and ki(t) are the phase space coordinates of the N test particles prop- 
agated with the Hamilton equations. As said, they are distributed like D{xi,pi,t). 
According to the theory of the Monte Carlo integration both integration procedures 
yield the same result in the limit of an infinite number of test particles. In practice 
one has to verify that the results do not depend on this number. Usually 100 test 
particles per physical nucleon in the system are considered as sufficient. It is very 
important to realize that these test particles have nothing to do with physical nucle- 
ons. They serve only as a representation of the one body Wigner density D{xi,pi, t). 
All observables which require more than its knowledge are beyond the scope of ap- 
plicability of these theories. Hence the possibility to extract source radii and hence 
correlation functions from the one body theories requires: 

• The smoothness assumption is valid 

• D{xi,Pi,X2,P2,t) ^ D{xi,pi,t)D{x2,P2,t) 

They are a consequence of the impossibihty to create two body Wigner densities or 
Wigner densities of two body observables like pi +p2 from the swarm of test particles 
defined as above. If both approximation were valid the correlation function is given 

by 
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— — = J d^xid^X2D{xi,pi)D{x2,P2){l ± cos{pi - P2){Xl - X2)) 



JlUp2 



= N(N-l) ^ COs{pi{to) - Pj{to)){Xi{to) - Xj{to))) (36) 

to is the (assumed common) freeze out time. The second approximation, the absence 
of two body correlations, is hard to control. The importance of many particle cor- 



relations for the fragment formation has been discussed in but its relevance for 
the proton or pion emission has not yet been investigated. 

VI. RESULTS FOR A GIVEN SOURCE DISTRIBUTION 

To interpret the different results given above it is useful to apply them to a situ- 
ation where the source is known. To keep the things simple we assume a completely 
chaotic source without any correlation between coordinate and momentum space: 

We start out from the Wigner density (eq.26) for a pair of particles as given in the 
QMD simulation. Averaging over the Gaussian source distribution we obtain for the 
correlation function as in eq.5 

C{p) = ^^TT^, = 1 ± (38) 

■> dpi dp2 ^ ' ' 

where S{k, r) = J S{ki, ri) ■ S{k2^ f2)d^Kd^R. If we apply the smoothness assumption 
(eq.6) we obtain 

r d^dwf±j^3. pg(k r)d^kd^r 

C^^iv) = "'^"^^^ ' = 1 ± e~P^(^+^) (39) 

f^^d^PS{k,r)dWr 

If we assume as in BUU, VUU or LV that the one body Wigner density is not given by 
Gaussians but as a sum over test particles (TP) each represented by a delta function 
in coordinate and momentum space 
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1 ^ 

where the p^s and distributed according to our source function we obtain 

as a correlation function 

C^^(p) = 1 ± e-^'^. (41) 

Defining the square of the source radii as j^^^f^g^ and comparing eqs. 38,39,41 
we observe that for the same measured correlation function C (p) we obtain different 
source radii depending on the simulation programs and the approximations used. 
From a mathematical point of view the difference between C'^^{p) and C^^{p) is easy 
to understand. Because the wave function used for the calculation of C^^{p) has a 
width of L, the true distribution of the source is the convolution of the distribution of 
the centers given by S{r, k) with the distribution of the one particle density around 
the centers. For C^^{p) one assumes that the true source distribution is given by 
S{f, k). In order to make both quantities comparable, both mean square radii have 
to be the same and hence one has to replace in C'^^{p) A by A'= L+A. 

However, being purely mathematical, this argument has an essential drawback. 
We have started out from the approximation (eq.l) that the source can be treated 
classically, and hence that the distance between two sources is large as compared to L 
0]. Hence, either the difference between A' and A is small and can be neglected. Then 
our approximation is valid. Or this difference is not negligible. Then our classical 
source approximation breaks down. That the wave function plays indeed a nontrivial 
role can be seen if one compares C^^{p) and C(p). In both cases the same single 
particle wave function has been employed. The result for the correlation function 
is, however, different. Only if L « A the difference between both is negligible. 
Hence our quantitative result confirms the well known qualitative argument that the 
smoothness assumption is only valid if the source can be assumed to be classical, i.e. if 
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the width of the wave function is small as compared to the size of the emitting system. 
Opposite, if L is of the same order as A as in all presently employed simulation models, 
the difference becomes important as we will see below and hence the smoothness 
assumption will break down. Hence we are confronted with the fact that present day 
simulation programs use a value of L which neither justifies the classical treatment 
of the source nor confirms the validity of the smoothness assumption. 

Nevertheless it seems that the community has agreed upon a pragmatic point of 
view in pretending that at least the classical treatment of the source is acceptable in 
modeling heavy ion collisions although a proof has not be given yet. Hence it may 
be useful to see whether under this assumption a quantitative prediction is possible. 
This includes the answer to two questions: To what precision we desire to measure 
the density and is the systematic error of the correlation function sufficiently small 
to obtain the desired precision. 

The study of the space time correlation is born out of the demand to measure 
the size of the system at the moment when the particles are emitted. If we study 
nucleons of the fireball, the density has to be in between twice and half normal 
nuclear matter density, because if the expanding fireball passes the latter density, 
there are no interactions anymore and hence the emission of particles defined as the 
time point of the last collision has ceased. For nucleons emitted from the spectator 
matter which remains at normal nuclear matter density one would like to know the 
source size. At lower energies the interest is to study whether the emitted nucleons 
come from a compound nucleus or whether the they are emitted from a subsystem 
called hot spot. Also here the density varies little around normal nuclear matter 
density. Whereas in the first case an uncertainty of the density determination of 
about 20% may be tolerable, the latter two require a precision of the determination 
if the source radius by about 3% (and hence of the mass number of about 10%) if one 
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would like to avoid that the uncertainty is already as large as the possible variation 
of the size of the system under investigation. 

In order to see whether this precision can be obtained we have to calculate the 
values for A, L and B for the cases of interest. If we assume that the rms radius 
of the source corresponds to the size of a nucleus at normal nuclear matter density 

1/3 

Rq = l-2Ajlf where Am is the atomic number of the nucleus we obtain 

A = A%'[fm'] (42) 

i.e. A = 21.5 fm^ for Am = 100 and A = 34 fm^ for Am = 200. Because our 
source emits particles according to a Maxwell Boltzmann distribution we can relate 
the slope in momentum space with the temperature of the emitting source and find 
that 

40 

, J /m'] (43) 

In the standard versions of QMD resp. IQMD the parameter L has the value 4.33 
and 8.66 /m^, respectively. Hence first of all we observe that L is not at all negligible 
as compared to A. However, as mentioned above, accepting a classical treatment of 
the source we can correct for this. It remains to be seen whether the smoothness 
assumption can be justified. Comparing the mean square radii obtained with and 
without the smoothness assumption 

RSA (L + A) 

= 1 hE. (44) 

we find that the smoothness assumption pretends a larger radius of the system. The 
value of F ranges between .87 for small systems at low temperature (5 MeV) and .98 
for large systems at high temperature (80 MeV). Hence for particles emitted from a 
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compound nucleus or from the spectator matter the error in the determination of the 
mass number due to the smoothness approximation is of the order of 20% even if the 
source is completely chaotic and of known form and the classical approximation of the 
source remains valid. For fireball nucleons the smoothness assumption produces an 
error of about 4% on the density. Of course if we were sure that we have a source of a 
given temperature we could also correct for the temperature, however such a source 
is not encountered in heavy ion physics where the excitation energy and hence the 
temperature changes in the course of time. 

VII. REALISTIC SIMULATIONS 

We have seen that for the most favourable condition (chaotic source of known 
form without any momentum space coordinate space correlation) the smoothness 
assumption enlarges the apparent source size by about 20%. If one applies now the 
simulation programs to real experiments one has to inspect the consequences of two 
facts: 

1) Nature most probably does not keep the rms radius of a nuclear wave function 
constant during a heavy ion reactions, QMD does. For observables which do not 
depend on the width of the wave function explicitly this may be of minor importance, 
the influence on observables which depend explicitly on the width, like the correlation 
function, is hard to judge since no calculations are available for a more sophisticated 
treatment of the reaction as done in QMD. If the width of the wave function has 
changed in the course of the reaction the difference between C^^{p) and C^^{p) 
cannot be corrected anymore by use of the known initial density distribution. 

In the QMD calculations the width of the wave function L serves two purposes. 
First it is used to have the proper one body density distribution when one initialize the 
nuclei. This is, however, a very weak condition because with much larger widths than 
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that actually employed one can obtain the same one particle density distribution. 
Second, it appears in the time evolution equations but only in form of the expectation 
value of the potential. Thus what counts for the time evolution is the convolution of 
the potential range and the width of the wave functions. Hence one can obtain the 
same expectation value of the potential for a smaller width and a larger potential 
range. Hence there is no need for an exact determination of the width L in the 
QMD calculation or, vice versa, the success of these calculations cannot be used to 
determine L. 

2) The source is as simulation programs show not at all chaotic and shows strong 
correlations between momenta and positions. 

Momentum space coordinate space correlations decrease the source size extracted 
from the correlation function as compared to the geometrical size of the source. 
This can be easily understood if one goes to the extreme. If the momentum is a 
monotonic function of the position, two particles with a small relative momentum 
have to come from the places very close in coordinates space. Thus the correlation 
function measures only that region in coordinate space from where these particles 
can come. Hence the stronger the momentum space coordinate space correlations are 
the smaller is the source size measured by the correlation function. As a consequence 
the value of A becomes smaller and the importance of the width of the wave function 
increases. Thus the stronger these correlations are the larger becomes the difference 
between R^^ and R (eq. 44). 

Thus for realistic calculations the situation becomes worse as compared to a static 
source. For a given size of the system correlations make A smaller and hence increase 
the importance of the width of the wave function if one compares eqs. 39 and 41 . 
They also do not allow to corrected for the smoothness assumption (eqs. 38 and 39) 
because the temperature is not anymore a global variable. The calculation of the 
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systematic error of the value of the radius determined by eqs. 39 or 41 requires more 
than the present models can predict, however the above arguments show that it will 
be larger than that for a static source. 

VIII. CONCLUSION 

We have discussed the possibility to extract source radii by comparing the ex- 
perimental results with the prediction of simulation programs. There is no doubt 
that the experimental results indeed show momentum space correlations caused by 
the bosonic or fermionic nature of the observed hadrons. These correlations carry 
information about the space time structure of the reaction. The goal to relate the 
observed correlation functions in momentum space with physical parameters in coor- 
dinate space like source radii, densities or energy densities can presently only achieved 
by use of transport theories. 

None of these transport models takes the bosonic or fermionic nature of the 
hadrons into account. Whereas this may be no essential drawback for many observ- 
ables it makes it impossible to calculate the correlation function in a straight forward 
manner. Every model requires the introduction of the (anti)symmetrization of the 
wave function in an ad hoc fashion in order to predict a correlation function. 

We have found that for all presently existing models, which can be subdivided 
into three classes, this introduction poses problems. 

Cascade models, in which classical particles are propagated, do not allow the 
calculation of a correlation function. The quoted values of source radii are totally 
artificial being a consequence of the employed approximation and not of physical 
origin. 

QMD, LV, BUU and VUU models allow the calculation of a correlation function. 
We find, however, that the basic approximation of the whole approach, namely that 
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the source can be considered as classical, is not fulfilled. Even if it were fulfilled, the 
systematic error of the extracted density introduced by the smoothness approxima- 
tion is for the most favourable case of a chaotic source of known form up to 20%. For 
realistic cases where space momentum space correlations are present and where we 
do not know the form of the source we have shown that the error will increase. This 
questions the possibility that in nuclear physics the HBT method allows a determi- 
nation of the density to a precision which allows to discriminate between different 
proposed interaction mechanisms. 

Of course this raises the question how to proceed. As we have seen we are plagued 
with systematic errors of the order L/A. There is first of all the open question whether 
the wave function of a emitted nucleon is smaller than L. Mean field calculation yield 
a much broader wave function and consequently the approximation of a classical 
source cannot be justified any more. Short range correlations, however, may distort 
this wavefunction. Hence it may be justified to address the question if there is a 
possibility to construct dynamical theories which can provide a prediction for the 
correlation function? Either one can try to decrease L or to avoid the systematic 
errors. 

The first suggestion implies a localization of the particles with a precision of about 
Ifm. This will be hardly possible. Not only because in a nuclear environment the root 
mean square radius of the wave function of the nucleon is considerably larger than 
the radius of a free nucleon but also because it implies an uncertainty of 200 MeV/c 
for the momentum of the nucleons which poses several severe technical problems for 
transport theories: 

• How to propagate particles in semiclassical theories whose velocity uncertainty 
is about 0.2c is unknown. 

• The sequence of collisions becomes undetermined 
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• The applied scattering cross sections have to be modified because the scattering 
partners are asymptotically not in a plane wave state. 

The second suggestion implies the construction of transport theories which prop- 
agates at least two particle wave functions and not parameters of the wave function. 
Presently such an approach is not available. 

Before a solution to these problems has been found the Hanbury Brown Twiss 
effect is a very nice quantal effect. Its application in nuclear physics to study the 
space time structure remains, however, premature. 

Interesting discussions with Drs. Ardouin, Erazmus, Gyulassy, Heinz, Lednicky 
and Werner are gratefully acknowledged. Furthermore I would like to thank Dr. 
Heinz for a careful reading of the manuscript. 
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